Electron correlation effects on the dielectric function of liquid metals 
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The acoustic excitations of the expanded metal solutions Li-NH3 have been measured by inelastic 
X-ray scattering as a function of the electron density by changing the Li concentration. The di- 
electric functions of these model metals with very low electron density have been derived from the 
high frequency sound velocity using the one component plasma approach corrected for screening. 
Their values, when combined with those from other metals, suggest that the electron correlation in- 
duced departure of the dielectric function from the random phase approximation follows a universal 
behavior whose main parameter is the electron density. 
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The study of the dynamical properties of interacting 
electrons is an important topic in solid state physics as 
shown by the great interest devoted to the understanding 
of phenomena as superconductivity and giant magneto- 
resistance. One way to obtain information on the dielec- 
tric response of a system of interacting electrons is to 
study the acoustic excitations of liquid metals, since the 
speed of sound in these systems is related to the screening 
action of the electrons on the ionic motion pj. In fact, in 
a simple Jellium model with the electrons screening the 
Coulomb forces between the ions, the speed of sound, c, 
can be expressed as: 



(i) 



Here Q is the wave number, e(Q) is the total dielectric 
function and fl pi is the ionic plasma frequency given by: 
^■pi — {^n^Ze) 2 ) /M, where M, Ze and rij are the mass, 
the charge and the number density of the ions, respec- 
tively. In this relation all the details of the electron dy- 
namics are hidden in the dielectric function, e(Q), that 
is purely electronic only if the ions contribution can be 
neglected - this is the case at high frequencies and Q- 
transfers, where the atomic response of the liquid is to- 
tally elastic and one measures the infinite-frequency value 
of the sound velocity, c. The strength of Eq. is to pro- 
vide a direct link between an electronic property as e(Q) 
and a ionic one as c. Therefore, the measurements of c 
offers the opportunity to check experimentally the valid- 
ity of theoretical approaches for e(Q) as, for example, the 
Random Phase Approximation (RPA), where 0: 
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Here, m is the electron mass, Dp the Fermi velocity and 
Qrpa is a measure of the inverse distance over which 
the self-consistent electric potential associated with the 
ions penetrates into the electron gas. This approxima- 
tion is exact in the limit of high electron density, or 
r s ^0. The parameter r s — [3/(47rn e )] 1 / 3 a ( ^ 1 , with ao 



the Bohr radius and n e the electron number density, is 
the ratio between the Coulomb and the kinetic energies 
of the electrons. The r s ^0 limit corresponds, therefore, 
to the weak coupling limit. As a matter of fact, the RPA 
works well in metals with a low value of r s and with no 
unfilled d or / shells, as shown in the case of Hg and Pb 
where r s is « 3 [2j. In metals with higher values of r s 
the RPA clearly fails as recently reported for the case of 
the Li-NHoSolution at the limit of solubility for which 
r s is « 7 0, 0, an d f° r heavy alkali metals It is 
believed that the failure of the RPA at large r s is mainly 
due to electronic correlation effects. In this context the 
Li-NH3 solutions provide an interesting model system to 
investigate the departure of the e(Q) from the RPA pre- 
dictions as a function of the electron density - here the 
r s value can be widely varied by changing the Li concen- 
tration from the saturated solution value r s —l A to the 
metal-insulator transition value ~ 12. 

In this Letter we report the experimental determina- 
tion of the high frequency velocity of sound in Li-NH3 so- 
lutions and the derivation of the corresponding dielectric 
function using Eq^ This provides a context to investi- 
gate the departure of the e(Q) from the RPA predictions 
as a function of the electron density - here the r s value 
has been changed between 7.4 and 10.4. These results, 
when complemented by literature data on liquid metals 
with lower values of r s 0, Q, H @, [H EJ H 111 13 El , 
provide i) the r st ^=4.6 value as the threshold value below 
which the RPA properly works and ii) the characteristic 
dependence of the electronic dielectric function on the 
electron density at r s larger than r s th- These findings 
suggest the remarkable result that the electron density is 
the relevant parameter in the description of the dielec- 
tric function even when the electron response starts to 
be affected by complex corrections to the RPA. 

The experiment was performed at the high resolution 
Inelastic X-rays Scattering (IXS) beamline ID16 at the 
European Synchrotron Radiation Facility in Grenoble 
(France). The measurements were carried out with a 
total energy resolution of ~ 1.5 mcV and an exchanged 
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momentum, Q, resolution of 0.4 nm -1 . Further details 
of the beamline are reported elsewhere ^(| . Every scan 
took about 210 min and the final spectra were obtained 
by averaging 3 to 4 scans. The Q values were selected 
between 1 and 15 ran -1 . The samples were prepared us- 
ing high purity Li (nominally 99.9 %) and high purity 
anhydrous ammonia (nominally 99.995%). A stainless 
steel container, with two sapphire windows (total thick- 
ness 500 /i) and with 5 cm 3 internal volume, was kept 
in a cryostat at T=240.0±0.2 K. The initial solution was 
prepared at the limit of solubility of lithium in ammonia 
(20±0.5 mpm). Successive dilutions were carried out in 
situ in order to measure three other metallic concentra- 
tions: 17.5±0.5 mpm, 9.3±0.2 mpm and 6.5±0.2 mpm. 

IXS spectra are reported in Fig. Q for the solution at 
17.5 mpm at a representative set of Q values. They con- 
sist of a central quasielastic line and two clear inelastic 
components which, as shown in Fig. disperse up to 
Q=7 nm -1 . The dependence of the measured spectra on 
the Li concentration is shown in Fig. at Q=2 nm -1 
for the four studied solutions. The spectral line shape 
changes considerably on changing the Li content: the 
quasi-elastic line becomes more and more intense with 
respect to the inelastic features as the metal to insula- 
tor transition is approached (top to bottom spectra in 
Fig. El). At larger Q values, particularly at the first sharp 
diffraction peak (Qm ~ 10 nnf 1 for the solution at 17.5 
mpm, see Fig. EJd), the spectra become narrow with a 
width comparable to the resolution function. 

The measured IXS spectra have been fitted to the con- 
volution of the resolution function with a model function 
composed of a Lorentzian curve for the central line and 
a Damped Harmonic Oscillator (DHO) function for the 
inelastic signal |l7j : 

qfn E/KbT f ! c (Q)T c (Q) 

Wl j l-e £ /^Hr2(Q) + £ 2 + 
I(Q)T(Q)n 2 (Q) x 
(n 2 (Q)-E 2 y+4T 2 (Q)E 2 )' 1 ' 

The quasielastic and inelastic spectral features are de- 
fined in terms of a set of two and three parameters, re- 
spectively. Those concerning the inelastic components 
are related to the energy, 7if2(Q), lifetime, r(Q), and 
strength, I(Q), of the acoustic-like excitations. Both 
Fig. and Fig. [5] show that the fitting function of 
Eq. describes well the experimental spectra. The 
Q-dispersions of the fitting parameter fi(Q) are shown in 
Fig. (left axes) together with X-ray measurements of 
the static structure factor, S(Q) 0] (right axes), and an 
indication of the location of twice the Fermi wavenum- 
ber, 2Kp (dotted vertical line). At low Q the acoustic 
excitations of the solutions show a linear dispersion (solid 
lines in Fig. as it is expected for a sound-like mode. 
In the panels (a) and (b) of Fig. El the results for the two 
more metallic solutions are reported. Here, f2(Q) shows 
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FIG. 1: Selected IXS spectra for the 17.5 mpm LiNH 3 solu- 
tion (o) at the indicated values of the exchanged wave vec- 
tor, plotted together with the total fitting function (-) cor- 
responding to Eq. |3] and the resolution convoluted inelastic 
features (dot-dashed line). At Q=l nm _1 and Q=10 nm _1 
the experimental resolution function (dotted line), scaled to 
the maximum of the spectrum, is also reported. 



a maximum between 5 and 6 nm -1 and then decreases 
reaching a minimum at a Q value close to the first maxi- 
mum, Qm, in the corresponding S(Q). This is a common 
finding in liquids which is usually referred to as the De 
Gennes narrowing |l7| . A different interpretation of this 
minimum formulated in terms of the Q-dependence of the 
electronic dielectric function has instead been proposed 
for the 20 mpm solution 0,0. In fact, for Q-values 
of the order of Kp, Eq. 0must be replaced by the more 
accurate Lindhard result Q which has a singular point at 
Q = 2 Kp. This singularity would in turn induce a kink 
in the dispersion curve via Eq. ^ However, given the 
proximity of 2Kp to Qm for both the 20 mpm and the 
17.5 mpm solutions, it is hard to disentangle structural 
and electronic effects. The dispersion curves of the solu- 
tions at 9.3 mpm and 6.5 mpm reported in the panels (c) 
and (d) of Fig.Elare, from this point of view, much more 
interesting since 2Kp and Qm are here considerably dif- 
ferent. As a matter of fact, these two dispersion curves 
both show a slight bending down in the proximity of the 
first diffraction peak (less evident for the solution at 6.5 
mpm) and a small kink at Q^8 nm , very close to 2 
Kp. This kink suggests indeed the presence of a Kohn- 
like instability of the system at Q ~ 2 Kp. This result - if 
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FIG. 2: Selected IXS spectra at fixed wave vector transfer 
Q=2 nm -1 for the four studied solutions. The experimental 
data (o) are plotted together with the experimental resolution 
(dotted line), the total fitting function (-) and the inelastic 
resolution-convoluted contributions (dot-dashed line). 



confirmed with a better statistics and higher Q-resolution 
- would then delineate the effects of the electron-phonon 
coupling in these expanded liquid metals, and support 
the interpretation proposed earlier 0, However, 
these data indicate that in any case the Kohn anomaly 
has only a small effect on the dispersion curve, and the 
minimum at Qm has mainly a structural origin. 

Linear fits to the dispersion curves at low Q, 
reported in Fig. are used to determine the 
high frequency sound velocity of the measured so- 
lutions: Cfixp (20 mpm)=1580±20 m/s, 0^(17.5 
mpm)=1519±60 m/s, CExp(9-3 mpm)=1430±60 m/s and 
CE£cp(6.5 mpm)=1450±60 m/s 20]. These values are 
compared with the RPA expectations obtained combin- 
ing Eqs. UlandU i.e. c RPA = v Fy /(Zm)/(3M) fl. In 
this expression, knowing that the Li ions form tetrahc- 
dral complexes with the ion surrounded by four ammo- 
nia molecules [^, |2^ , M is identified as the mass of the 
Li(NH 3 )4 complex, as already proposed in refs. pi HfSj. 
In this picture, the solution at the limit of solubility (20 
mpm) is completely made up of these complexes, while 
at lower Li content, the solutions can be viewed as a 
two component system where the Li(NH3)4 complexes 
are still stable entities embedded in the uniform back- 
ground of the ammonia molecules [23j | . The comparison 
between the CExp and crpa is presented in Fig. 4 as 
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FIG. 3: Dispersion curves for the acoustic excitations ob- 
tained from the fit procedure discussed in the text for the 
Li-NH3 solutions at the indicated concentrations. The X-ray 
S(Q) spectra are also reported (right axes) |T3 |. The vertical 
dotted lines indicate the wave vector position corresponding 
to 2 Kf- The solid lines are the low Q linear fits to the 
dispersion curves. 



also report the results obtained for other liquid metals, 
as taken in the literature from ref. for Al, ref. 
for Hg, ref. H3 for Pb, ref. [HI for Li, ref. 07J for Na, 
ref. Q for K, ref. [Hf for Cs and ref. @ for Rb. At 
low r s , as expected, a very good agreement exists be- 
tween the experimental results and the RPA predictions, 
confirming nicely that Eq. prelates the infinite frequency 
sound velocity and the dielectric function. A clear depar- 
ture from the RPA is observed in Fig. 4 at the threshold 
value r s t/i=4.6, above which CExp becomes sensibly larger 
than crpa- 

The very intriguing result of Fig. 4 is that all the ex- 
isting data follow a continuous trend despite of the very 
large variety of considered systems: steric effects, chemi- 
cal specificities, and spatial extent of the different atoms 
and complexes are somehow irrelevant for c Exp / c RPA 
and eppAl '(-Exp - the only important parameter is r s . In 
other words, regardless of a number of parameters not 
correctly treated by the RPA, as electronic correlations 
and finite ionic size effects, the electron dynamics has a 
general dependence on the electron density. 

The results reported in Fig. 4 are further emphasized 
considering the r s -dependence of the screening wave vec- 
tor QExp—^pi / CExp, shown in the insert of Fig.0] Here, 
Qexp is compared with Qrpa and with the SSTL theo- 
retical prediction |24| , which is an attempt to extend the 
range of the RPA by introducing a local field correction 
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FIG. 4: Ratio between the experimental and the RPA squared 
sound velocities (left axis) as a function of r s for different liq- 
uid metals 0, IS E3, E3, E3, E3, E3, E3l These data correspond 
via Eq.l to the ratio between the RPA and the experimental 
dielectric function (right axis). The dotted line corresponds 
to the unity value. In the insert the corresponding experi- 
mental screening wave vectors (closed circles) are compared 
to the RPA (solid line) and to the SSLT theory (dot dashed 
line). The dotted line emphasizes the phenomenological rj 3 
law obeyed by metals for r s >4.6. 



readjusted to an external field. Once more, experiment 
and theories are in good agreement up to r st h, while a 
clear deviation is observable above r st h- We also note 
that the SSLT theory not only does not improve RPA in 
describing Qexp, but seems to predict a departure from 
RPA opposite to the one experimentally measured. Fi- 
nally, it is very interesting to note that for r s > r st h, 
Q Exp is extremely well reproduced by an empirical r s 3 ^ 2 
dependence. 

In conclusion, we presented IXS measurements in ex- 
panded liquid metals with high r s allowing the determi- 
nation of the high frequency sound velocity: these data, 
along with literature data for metals with lower r s , em- 
phasize the breakdown of existing electron gas theories 
above a threshold value r st h- Most interestingly, the di- 
electric function, when the electron-electron correlations 
become increasingly important, shows a departure from 
RPA that follows a characteristic behavior which depends 
only on the electron density and not on the specific de- 
tails of the considered systems. Future theories on the 
electronic dynamics must be compared with such a find- 
ing in order to test their ability to describe the deviations 
from the RPA. 
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